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Abstract 

An hallmark of present-day Dyson-Schwinger equation applications in hadron 
physics is the existence of a systematic and symmetry preserving truncation 
^ ' scheme. This enables the proof of exact results; e.g., the leptonic decay con- 

^ . stant of every pseudoscalar meson except the pion vanishes in the chiral limit. 

\ Calculations using the scheme's leading-order truncation are reliable in the 

^ ■ vector and flavour nonsinglet pseudoscalar channels. In this rainbow-ladder 

O . truncation, an impulse approximation provides the consistent current for all 

^ \ six-point quark Schwinger functions. That is well illustrated via the anoma- 

O ' lous process vr'^ — > 77. Using two-, three- and four-point Schwinger func- 

^ ■ tions calculated in rainbow-ladder truncation, the textbook value of the width 

^ ! is obtained algebraically and independent of model details if, and only if, the 

' impulse approximation is used to describe the associated matrix element. 

cn : 

O 1. INTRODUCTION 

^ \ A number of brief accounts of the modem application of Dyson-Schwinger equations (DSEs) to hadron 

■ physics and QCD have recently appeared 1 1 2»i3j|4,,5|. They have recapitulated on aspects of dynam- 

^ . ical chiral symmetry breaking and confinement, and the expression of these phenomena in the hadron 

'-^ \ spectrum; and contain original contributions regarding: comparisons between DSE results and lattice 

^ ■ simulations, and the mutual benefits of such comparisons 1 1 4|; the pion and its excited states |2||3l; and 

^ . states with vacuum quantum numbers |5 1. We have a modicum to add. 

c : 

> ; 2. GAP AND BETHE-SALPETER EQUATIONS 

^ . We begin with QCD's renormalised gap equation 

' A a 

- ' S{p)-^ = Z2{ij-p + mbarc) + ZiJ g^D^,{p - q)^^^S{q)mq,p) , (1) 



wherein: D^y{k) is the dressed-gluon propagator; T1{q,p) is the dressed-quark-gluon vertex; m-bare is 
the A-dependent current-quark bare mass; and jj^ := d'^q/ (27r)^ represents a translationally-invariant 
regularisation of the integral, with A the regularisation mass-scale. We employ a Pauli-Villars scheme, 
which is implemented in Eq. Q by considering the quarks as minimally anticoupled {g^^ = ig) to 
additional massive gluons {m^^ = A). This effects a tempering of the integrand that is expressed via a 
modification of the gluon propagator's ultraviolet behaviour: 

1 1 1 

(2) 



{p — q)^ {p — q)^ {p — qY + ' 

through which the integral's superficial linear divergence is regulated. The regularising mass-scale is 
naturally removed to infinity as the last stage in the calculation of any renormalisable quantity. In 
Eq. O, Zi 2(C^7A^) are respectively the quark-gluon-vertex and quark wave function renormalisation 
constants, which depend on A and the renormalisation point, C„ as does the mass renormalisation con- 
stant ZmiC^?) = ^4(C^, A^)/2'2(C^, A^). The gap equation's solution has the form 

S{p)-^ = i7 .pA{p\ e) + B{p\ e) ^ YW^^ [il-p-r M(/)] , (3) 



where M{(^) = m{Q) := mbare(A) Z^^(C^, A^) is the running quark mass,' and is obtained subject to 
the renormalisation condition 



S{p) 



-1 



■ p + m{() 



(4) 



The renormahsation scheme should be understood to yield mass-independent renormalisation constants. 

The gap equation illustrates the features and flaws of each DSE. It is a nonlinear integral equation 
for S{p) and hence can yield much-needed nonperturbative information. However, the kernel involves 
the two-point function D^i,{k) and the three-point function T'^{q,p). The gap equation is therefore 
coupled to the DSEs these functions satisfy. Those equations in turn involve higher n-point functions 
and hence the DSEs are a tower of coupled integral equations with a tractable problem obtained only once 
a truncation scheme is specified. It is unsurprising that the best known truncation scheme is the weak 
coupling expansion, which reproduces every diagram in perturbation theory. This scheme is systematic 
and valuable in the analysis of large momentum transfer phenomena because QCD is asymptotically free 
but it precludes any possibility of obtaining nonperturbative information, which we identified as a key 
aspect of the DSEs. 

The gap equation describes the manner in which the propagation of a quark is affected by the 
medium being traversed, and the same is true of the gauge sector gap equations satisfied by gluon and 
ghost Schwinger functions. The Bethe-Salpeter equations are the necessary complement as they provide 
the means by which one may understand how QCD's elementary excitations form the colour singlet 
bound states that are the observable manifestation of strong interaction dynamics. 

An important example is the inhomogeneous Bethe-Salpeter equation for the axial-vector vertex: 



{k-P) 



tu 



+ 



- tu 



(5) 



where X5fi{q; P) = S{qi)T5^{q; P)S{q2) with S{£) = diag[S'„(£), 5rf(£)] and, plainly, we focus on 
SU{Nf = 2) so that {T^;j = 1,2,3} are Pauli matrices in flavour space. The conventions of Eq. (|5ll 
describe a vertex with total momentum P and relative momentum k, whose outgoing quark and antiquark 
legs carry momenta ki = k + rjP and —k2 = —k + {l — r])P, respectively. The appearance here ofr]£ 
[0, 1] expresses the arbitrariness in defining the relative momentum in a Poincare covariant framework: 
physical observables will be independent of ry in a symmetry preserving truncation, as seen, e.g., in 
Ref. E|. Of particular importance in Eq. ^ is K{q, k; P), the fully amputated dressed quark-antiquark 
scattering kernel. This four-point Schwinger function also appears implicitly in QCD's gap equation 
because it is one piece of the kernel in the inhomogeneous integral equation satisfied by r^{q,p). 

3. PION AND RELATED MATTERS 

A primary measure of our comprehension of strong QCD is the possession of a veracious understanding 
of the pion, and for this a clear grasp of the nature of chiral symmetry and its dynamical breaking are 
vital. In asymptotically free theories, chiral symmetry and its breaking are completely expressed via the 
axial-vector Ward-Takahashi identity: 

P^ri^{k;P)=S-\k+)iro^ + ij,^S-Hk^)-iM{C)ri{k;P)-Ti{k;P) (6) 

wherein the new elements are the current-quark mass matrix: A4{C) = diag[mu(C), ^rid the 
pseudoscalar vertex 



m;P) 



tu 



75y 



tu 



+ / [X6iq;P)lrKi:{q,k;P). 



(V) 



'The dressed-quark mass function is renormalisation-point-independent; i.e., M{p^ , — M{p'^ , ("^ 



The last two terms of Eq. Q vanish in the chiral limit. The identity is then a concise statement of 
chiral symmetry, a consequence of which is plainly an intimate relation between the three-point function 
in Eq. ^ and the two-point function in Eq. Q- Since dynamical chiral symmetry breaking is an essen- 
tially nonperturbative phenomenon, satisfying our measure of comprehension requires that there exist a 
nonperturbative and systematic truncation scheme which preserves identities such as Eq. Q. 



3.1 Model-independent results 

At least one such scheme exists |71|8|: a dressed-loop expansion of the dressed-quark-gluon vertices that 
appear in the half-amputated dressed-quark-antiquark scattering kernel: {SS)K, and that is sufficient to 
arrive via the DSEs at exact results in QCD. For example, it is a general property that the axial-vector 
vertex, obtained via Eq. (|5l), exhibits a pole whenever = —m^^, where mjr^ is the mass of any 
pseudoscalar u-d meson; viz.. 



Ti(k;P) = regular terms + J^'^^^^ (k;P), (8) 

" n 7Vn 



with rl.^{k; P) being the + bound state's Bethe-Salpeter amplitude: 

Ti^ik; P) = r^75 [iE^„ {k; P) + 7 • PF^^{k; P) + j . k k ■ P G^„{k; P) + k^P, H^„{k; P)] , 

(9) 

and /tt^ , its leptonic decay constant 

U„6'^P^ = Z2tr riTS7/.5(g+)rijg;P)5(g_), (10) 

Jg 

where the trace is over colour, flavour and spinor indices. Equation ( fTUb is the precise expression in 
quantum field theory for the pseudovector projection of the meson's Bethe-Salpeter wave function onto 
the origin in configuration space. The appearance of Z2 guarantees that /,r„ is gauge invariant, and 
independent of the regularisation scale and renormalisation point. 

The lowest mass pole contribution to Eq. denoted by n = 0, is the ground state pion, which 
receives the spectroscopic assignment N "^^^^L j = 1 ^Sq in the naive qq quark model. The n > 1 pseu- 
doscalar meson poles correspond to those excited states of the pion that would receive the spectroscopic 
assignments (n + 1) in the quark model. 

The pseudoscalar vertex, T\{k; P), which appears in Eq. also exhibits such a pole: 

iTl{k- P) = regular terms + ^ , rij/c; P) , (11) 

and here the residue 

ip^^iOS'^ = Z,tT rirS5(<z+)rijg;P)5(g_) (12) 

Jg 

expresses the pseudo^caZar projection of the meson's Bethe-Salpeter wave function onto the origin in 
configuration space. The appearance of Z4 guarantees that p^,^ is gauge invariant, independent of the 
regularisation scale, and evolves correctly under the renormalisation group. 

In QCD it therefore follows from Eqs. ^ and (fTTT l. as a necessary consequence of chiral symmetry 
and its dynamical breaking, that for any pseudoscalar u-d meson t2i|3j|9J 

U^ml^ = K(C) + m,(C)] ■ (13) 

(This expression is valid independent of the magnitude of the current-quark mass and, following Ref. |6l, 
the generalisation to SU{Nf) -flavour is straightforward.) The so-called Gell-Mann-Oakes-Renner rela- 
tion for the ground state pion appears as a corollary of Eq. (IT3l llUnOl, and another important corollary 



of Eq. (fT3l . valid for pseudoscalar mesons containing at least one heavy-quark, is described in Refs. 1111 . 
In addition, the proof of Eq. (fT3b establishes |2 1 that in the chiral limit 

/° ^„ := lim , = 0; (14) 

namely, in the chiral limit the leptonic decay constant vanishes for every pseudoscalar meson except the 
ground state pion. This result in particular is a constraint on models and theoretical methods used to 
search for exotics and hybrids. 



3.2 Ah initio calculations 

The leading-order term in the scheme of Refs. |2l IHl is the renormalisation-group-improved rainbow- 
ladder truncation, which is expressed in the Bethe-Salpeter equations via 



Ki:iq,k;P) = -47raiQ')D'^:%Q) 



Y7p 



ts 



y7. 



(15) 



with a{Q'^) an operative coupling, and in the truncated gap equation through: 

S{p)-' =Z2{h-p + mbarc) + 47ra((p - qf) D'^l%p - q)^j^S{q)^j, . (16) 

This truncation can be used to correlate lattice data on dressed-quark and -gluon Schwinger func- 
tions fT2l . It is also a reliable tool in the study of vector and flavour nonsinglet pseudoscalar channels 
1 8 1, and has been refined and exploited by Maris and Tandy in a series of articles^ via a one-parameter 
model for the operative coupling: 



Q2 CO' 



r+(l + QVA2jCD 



y^(Q^), (17) 



wherein JP'(Q2) = [i _ exp(-QV nit = 0.5 GeV; r = - 1; 7^ = 12/25; and Aqcd = 

0.234 GeV. This form expresses the interaction as a sum: the second term ensures that perturbative 
behaviour is preserved at short-range; and the first makes provision for enhancement (strong coupling) at 
long-range. The true parameters in Eq. (fTTt are D and uj, which together determine the integrated infrared 
strength of the rainbow-ladder kernel. However, they are not independent: in fitting to a selection of 
observables, a change in one is compensated by altering the other; e.g., on the domain uj € [0.3, 0.5] GeV, 
the fitted observables are approximately constant along the trajectory lo D = (0.72 GeV)^ fTSl . Hence 
Eq. (fTTb is a one-parameter model. This correlation: a reduction in D compensating an increase in uj, 
ensures a fixed value of the interaction's integrated infrared strength. It also defines a single dressed-glue 
mass-scale that characterises infrai"ed gluodynamics: 

mg = 720MeV. (18) 

Using Maris and Tandy's model for the operative coupling, the rainbow-ladder truncation yields 
fBl /tto = 0.092 GeV; m^^ = 0.14 GeV; p^^, = (0.81GeV)2, at a current-quark mass md(lGeV) = 
mu(l GeV) = 5.5 Me V. This model is currently being employed in a rainbow-ladder analysis of n 7^ 
pseudoscalar mesons, and it has been reported |2j|3 that when only the ET^^{k; P) term in Eq. ^ is 
retained,^ one obtains: m^, = 1.1 GeV cf. m°^P*- = 1.3 ± 0.1 GeV; and /^^^^ < 1.5 MeV at the 
physical (nonzero) current-quark mass. The latter result indicates that Eq. (fT4t is a good approximation 
even at real-world current-quark masses and thereby emphasises its usefulness as a constraint. 

^An exemplary success was their prediction 1 13 1 of the electromagnetic pion form factor 1141 . In addition, see Refs. (BlIlH 
and references therein, and the review of their substantial contributions in Ref. 1171 . 

''With this amplitude alone one obtains values for m^p, /^(, that are approximately 25% smaller than experiment. 



4. ABELIAN ANOMALY 



A correct description of the neutral pion's two-photon decay is another stringent test of one's under- 
standing of strong QCD, and that understanding again relies fundamentally on a true comprehension of 
chiral symmetry and its dynamical breakdown. To illustrate this we consider the renormalised impulse 
approximation to the axial-vector-vector-vector (AW) vertex: 

eM 

T^^,p{ki,k2)=iT 5(4+)ri^(4+,^-o)5(^_o)iQr^(^-o,^)5(^)iQr,(^,£o+), (19) 

where l^p = ^ + a^i + /9^2 and the electric charge matrix Q = diag[eu, e^] = ediag[2/3, —1/3]. In 
Eq. ( fT9l . the dressed-quark propagator is determined by the rainbow gap equation, and the dressed vector 
and axial-vector vertices are obtained from the dressed-ladder truncation of the relevant inhomogeneous 
Bethe-Salpeter equations. It is apparent from Eq. ([Sll that for {ki + ^2)^ + k, 0; namely, in the 
neighbourhood of the pion pole, T'^yp{ki , ^2) is dominated by the vr*^ 77 coupling and therefore plays 
a role in describing that decay. 

The bare AW vertex exhibits a superficial linear divergence and as with all other Schwinger 
functions it must be rigorously defined via a translationally invariant regularisation scheme. In this case 
the appropriate Pauli-Villars prescription corresponds to minimally anticoupling the photon to additional 
flavoured quarks with a large mass m^^ = M. To elucidate we introduce 



T'^,p{ki,k2;rn) := tr 5^(4+) ^^"(4+, ^-0) ^^(^-o) iQr™(£_o, ^) iQr^Tl^, , (20) 

where m is the renormalisation point invariant current-quark mass and every element in the integrand is 
calculated using this cunent-quark mass in the appropriate rainbow-ladder truncated DSE. In this case 
Eq. ( fT9t is explicitly 

r3^^(/ci, /C2; m) = fl,p{kr, k2;m) - f^^pik^M] M) , (21) 

with M ^ 00 as the last step in the calculation. 

We now focus on the chiral limit, m = 0, and consider the divergence of the AW vertex at 
= {ki + A;2)^ ~ 0; viz., in the neighbourhood of the massless pion's pole: 



PpT^,p{kiM; 0) p,^^^^^^^^,^^ = PpT'^,p{kiM\ ^) - PpT^,p{kiM: M) (22) 
= Z2(C^ tr ^757.Qy { 5^(^-o) iQr™(£_o,£) " 5m(^-o) i QT;!^ (^-0 , ^) Sm{1)} 

- Z2(C^M2)tr^757;.Qy {5AWiQC(£,£o+)5A(4+) -5A/WiQr*^(^,4+)5A/(4+)} 

+ 2M(C)tr^5Af(4+)^ri*^(4+,^-o)5M(^-o)*Qr;i'^(^-o,^)5MWiQr^(^,4+) • (23) 

Equation (123 1 follows from the vector vertex analogue of Eq. ^ and the axial- vector Ward-Takahashi 
identity depicted in Fig.fn"* which is an extension of 

PpS{k^)Ti^{k-,P)S{k.) = i^J-^S{k.)+S{k+)i^J-^-S{k+){M{0,^ (24) 

and can be derived following the method in Ref . ITOl : the identity is valid if, and only if, every dressed- 
quark propagator that appears therein is obtained from the rainbow DSE, Eq. (fT6b . and the accompanying 
dressed vertices are determined from the consistent ladder Bethe-Salpeter equation; i.e., the equation 
constructed using Eq. (I15t and the rainbow quark propagator. It is plain that in the chiral limit Eq. d23l 
does not contain terms which depend on the light-quark mass. 

*The relevance of this Ward-Takahashi identity is plain if one re-expresses the vector vertices through: = 
^27ji G (SS)^^, where G is the unamputated renormalised rainbow-ladder-dressed quark-antiquark scattering matrix. 
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Fig. 1: This axial-vector Ward-Takahashi identity is an extension of Eq. <24> . It is valid if, and only if: the dressed-quark 
propagator, S, is obtained from Eg. <16> : the axial- vector vertex, Fs^, is obtained from Eq. Jsj with the kernel constructed 
from S and Eq. <15t : the pseudoscalar vertex is constructed analogously; and the unamputated renormalised quark-antiquark 
scattering matrix: G = (SS) + {SS)K{SS) + {SS)K{SS)K{SS) + [. . .], is constructed from the elements just described. 



Consider the first term on tiie r.li.s. of Eq. (I23t . The compiete poiarisation integral is convergent, 
by virtue of our regularisation, and must yield a rank-two Lorentz pseudotensor depending on only one 
four- vector argument, ki . A nonzero value is impossible. The second term vanishes for the same reason. 

Let's turn to the last term. The only relevant contribution appears in the limit M oo, in which 
case all dressing is suppressed and we ai^e left to evaluate 



T^^{ki,k2;M) 



. d'^l I I 1 

/ TTru -■ — -ff — n^7 V -■ — ff — ^2^^ -■ — T-irf ^Q^^ (25) 

1 1 1 



f d £ 111 
{el - el) M tvD / — rj ^75 ^ — rj J—TJ ^ly (26) 

^M-P- f'lP / ^ (£2^ + M2) (£2 ^ + M2) (^2 + ^2) ^^^^ 



= i- e^upakipk2a / (ix / dy—^ — -^-^ • (28) 

vr Jo Jo M2 + xy (fei + fc2)^ 

Recall now that at any order in the nonperturbative, systematic, symmetry-preserving DSE trun- 
cation scheme of Refs. IH, dynamical chiral symmetry breaking entails that in the chiral limit the 
divergence of the axial-vector vertex is dominated by the pion pole in the neighbourhood of = 0. 
This result is expressed in Eq. (ISjl. It follows that in the rainbow-ladder truncation, which is the leading- 
order in that scheme. 



p2 = {A;i+fc2)2faO 



Je 

Consequently, in the chiral limit, the Bose-symmetrised amplitude for vr*^ —>■ 77 decay 

Q; 
jTT^^-^-rr = T^u{ki,k2) + r,,^(A;2, ki) = i—- £p,^p^ kip k2a ■ (30) 

In Eq. d30l we have arrived at the textbook result with no dependence on model details; viz., in- 
dependent of the detailed form of the operative coupling, Eq. (fTTl . and thereby provided an explanation 
of the numerical analysis reported in Ref. |20|. In doing so we have made plain that the impulse ap- 
proximation expressed in Eqs. ( fT9t and ( l29l is the precise and only truncation of the six-point Schwinger 
functions relevant to this decay which is consistent with the rainbow-ladder truncation. It follows that the 
"seagull terms" identified in Ref. 1211 are absent at this order but may appear in a consistent treatment of 
the process at an higher order of truncation. Moreover, it is generally true that seagull terms are absent 
in the consistent rainbow-ladder analysis of any six-point quark Schwinger function. 



5. EPILOGUE 

There are many additional applications of interest to this audience, of which a variety are reviewed in 
Refs. 1 17 22 1, and one may say that the keystones of success in modern DSE usage are: an appreciation 
and expression of the momentum-dependence of dressed-parton propagators at infrared length-scales; 
and the identification and employment of a nonperturbative, systematic, symmetry preserving truncation 
scheme. 

One topic of material importance, worth mentioning explicitly here, is the pion's valence quark 
distribution. The DSE calculation 1231 agrees with perturbative QCD's prediction |24|; i.e., nj^(x) cx 
(1 — x)^ in the valence region,^ but this is inconsistent with extant data l25l and a recent transverse lattice 
calculation |26|. The discrepancy between data and the expectation from perturbative QCD cannot be 
resolved using contemporary simulations of lattice-regularised QCD 127J. It raises difficult questions 
and has prompted proposals for remeasurement of the distribution function l28l . 

We have omitted a discussion of contemporary challenges confronting the application of DSEs 
in hadron physics and QCD. Important amongst these are: the development of a veracious picture of 
the infrared behaviour of the Bethe-Salpeter kernel; and the extension of the framework to processes 
involving baryons. This list is not exhaustive. Nonetheless, comments on the former may be found in 
Refs. Hinil^llOlEIl and the latter in Refs. 1321 . 
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